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I. INTRODUCTION 



Many of the theoretical and phenomenological analyses of QCD processes rely upon the 
factorization of the underlying dynamics into a short-distance dominated part, amenable to 
QCD perturbation theory, and a large-distance part that has to be taken from experiment 
or be determined by nonperturbative methods. Among such processes, the pion form factors 
(electromagnetic and transition) play the role of a theoretical laboratory to test various ideas 
and techniques. The key ingredient in these descriptions is the pion distribution amplitude 
(DA) ipT,{x) which represents the pion bound state. At leading twist two it is defined in 
terms of a nonlocal axial current and reads [1] 



(0|(i»7^75[^,0]M(0)|7r(P))|, 



= ^/.P^ /'rfxe^^(^-^V^'n^,/io), (1-1) 

^0 



where x (x = 1 — x) is the longitudinal momentum fraction carried by the valence quark 
(antiquark) in the pion and the path-ordered exponential, i.e., the lightlike gauge link 



[2,0] = Pexp 







;i.2) 



ensures gauge mvariance. 

The pion DA has an expansion in the basis of the Gegenbauer polynomials which con- 
stitute the eigenfunctions of the one-loop meson evolution equation. At a typical hadronic 
scale /Iq, which serves as a normalization scale, one then has [2] 



V9(*=2)(x;/i^) = V.-(x) 



1 + a2(/io) Cl'\2x - 1) + a4(/io) '(2x - 1) + . . . , (1.3) 



where (p^{x) = 6xx is the asymptotic pion DA. By virtue of the leptonic decay vr — )■ jJi^ffj,, 
one obtains the normalization dx (pn~'^\x, jj^) = 1, which fixes Oq = 1. 

Rather than try to derive the pion DA as a whole, one attempts to reconstruct it from 
its first few moments 





where ^ = 2x — 1. The values of the moments may be determined by means of QCD sum rules 
(SR)s with local [3] or nonlocal condensates [4-7], or be computed by numerical simulations 
on the lattice [8-10]. Once they are known, one can use them to reverse engineer the pion DA, 
with a precision depending upon the influence of the magnitude of the discarded higher-order 
moments. It was shown in Ref. [11] that, using QCD sum rules with nonlocal condensates, 
one can de facto resort to the first two Gegenbauer coefficients 02, 04, while aj with i = 6, 8, 10 
turn out to be negligible. Once the shape of the pion DA has been determined at some (low) 
normalization scale /i^ around 1 GeV^, one can evolve the Gegenbauer coefficients to higher 
values of the momentum scale using the Efremov-Radyushkin-Brodsky-Lepage [2] evolution 
equation which is determined by means of QCD perturbation theory. 

It turns out that another quantity which is intertwined with the form factors of the pion 
is its inverse moment 

/•I 1 

(x = dx - ip^{x) . (1.5) 
Jo X 

This quantity is one of the key ingredients of the pion-photon transition form factor, a process 
that has attracted the continuous attention of theorists [11-27] and experimentalists [28- 
30]. Actually, the most recent measurement of this observable by the BaBar Collaboration 
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[30] has provided controversial results. At moderate values of the momentum transfer, up 
to 10 GeV^, the new high-precision BaBar data agree well with the previous CLEO data 
[29] and can be best described by pion DAs that have their endpoints strongly suppressed 
[27, 31, 32] — as the Bakulev-Mikhailov-Stefanis (BMS) model [11] derived from QCD sum 
rules with nonlocal condensates. By contrast, the high-Q^ BaBar data show an unexpected 
growth with which cannot be understood on the basis of collinear factorization and calls 
for pion DAs that have their endpoints strongly enhanced [33, 34]. This intriguing behavior 
consists the basic motivation for the present investigation, though we will not attempt to 
describe any data. 

We shall employ in this work the method of QCD SRs with nonlocal condensates (NLC)s 
with the aim to estimate the slope of the pion DA in the region a; ~ 0, trying to understand 
the fine structure of the pion DA in this region vs. the ansatz for the quark-virtuality 
distribution in the nonperturbative QCD vacuum. Our main interest will be in the behavior 
of pion DAs with distinct endpoint characteristics. QCD SRs were mainly proposed with the 
purpose of studying the integral characteristics of the pion DA. To overcome this restriction, 
we shall design an operator for defining integral derivatives of the pion DA. These will 
supplement the results obtained with SRs which employ the standard derivative of the pion 
DA. In this latter case, we will use in our analysis not only a delta-function ansatz for the 
vacuum quark-virtuality distribution, but also a refined model which describes the large- 
distance regime more accurately. 

The paper is organized as follows. In Sec. II, we briefly discuss the main features of the 
nonperturbative QCD vacuum in terms of the scalar quark condensate. We also give the 
corresponding SRs and recall their main ingredients. Section III is devoted to the calculation 
of the slope of the pion DA in the endpoint region employing two different techniques: integral 
SRs and differential SRs. Finally, Sec. IV contains our conclusions, while some important 
technical details are given in four appendices. 

II. NONPERTURBATIVE QCD VACUUM WITH NONLOCAL CONDENSATES 

The basic idea underlying the NLC approach is that the vacuum condensates possess 
a correlation length which endows the vacuum quarks with a non-zero average virtuality 
(kg) (see, for instance, [35]). To analyze the nonlocality of the vacuum condensate, it is 
useful to parameterize the lowest one* (g(0)[0, z]q{z)) = Ms{z'^) with the help of the vacuum 
distribution function /^(a): 

oo 

Ms{z') = (qq) j f s{a) e'^^' ' Ua (2.1) 



that describes the distribution of the vacuum-quark virtuality a [4]. Assuming fs{ot) = 
6{a — \q/2), that takes into account only a fixed virtuahty of the vacuum quarks, leads 
to the simplest Gaussian model 

(g(O)g(^)) = (g(0)g(0))e-l^^l^?/« (2.2) 

for the scalar quark condensate [4]. 



* In this work we use the gauge z'^A^ = 0. Therefore, one has [0, z] = 1. 
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The parameter represents the typical quark momentum in the vacuum and is given by 



mm) - ^- 

In this work, we use the value A^ = 0.4 GeV^, which is supported by several analyses, 
though values within the interval [0.35-^0.45] GeV^ are still acceptable (see [11, 22, 36] and 
references cited therein). 

The QCD SRs with nonlocal condensates for the pion DA were first proposed in [4] and 
were significantly improved in [11] from which we quote 



oo 

o2 /n#2 



so 

+A$G(a;, M^) + [A$s(a;, M^) + A$v(2:, M^) + A$T(a:, M^)] ^ . (2.4) 

Here, ipAi is the Ai-meson DA, whereas f^, and are, respectively, the decay constants of 
the Ai and the vr-meson. The Ai-meson state is an effective state that collects the ir' and 
the ai meson. The nonperturbative ingredients in the theoretical part of the SR are the 
gluon-condensate term A$g'(x, M^) and the quark-condensate contribution [...]q. This lat- 
ter contribution contains the vector-condensate term (V), the mixed quark-gluon condensate 
term (T), and the scalar condensate term (S). The explicit expressions for the nonperturba- 
tive contributions and the NLO spectral density p^^lf^ (x) are given in Appendices A and 
C, respectively. It turns out that in the endpoint region, the first radiative correction in 
the spectral density, which is of 0{as), is too large, thus overshadowing the zeroth order 
perturbative contribution and the nonperturbative contribution. For that reason, we use in 
this work the leading-order (LO) approximation p^peft{x) = 3xx/27r'^. In order to include 
radiative corrections into the spectral density — when analyzing the endpoint region — one 
would have to resum all radiative corrections, a formidable task outside the scope of the 
present investigation. 



III. SLOPE OF THE PION DA 

As already mentioned in the Introduction, the endpoint region of the pion DA turns out 
to be of particular importance for a variety of pion observables — see [16, 37, 38] for an in- 
depth discussion of this issue. For example, in the case of the pion-photon transition form 
factor, one finds that in order to comply with the CLEO data — and those BaBar data close 
to them — one needs endpoint-suppressed pion DAs, while the high-Q^ BaBar data can only 
be described with flat-type pion DAs. Therefore, it is crucial to explore the fine details of 
the pion DA in a region around the origin a; ~ 0. One way to study the endpoint regime of 
the pion DA is provided by the inverse moment [39]. In Fig. 1 we show the values of 

{x-\y)) = \fj'^-^dx (3.1) 

with 5 = 0.05 for some pion DA models with a characteristic behavior in the endpoint region. 
From this figure one sees that the BMS pion DA (dashed green line), derived with the aid of 
nonlocal condensates [11], exhibits an evident endpoint suppression, while all other models 
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have more (flat-top — dashed-dotted-dotted red line — and CZ — dashed-dotted green line) or 
less (asymptotic DA) endpoint enhancement. The considered models of the pion DA, and 
their parameters, are given in Appendix B. In the present work we will probe the endpoint 
region of the pion DA also by another means, namely, the "integral" derivative which will 
be defined next. 
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FIG. 1: Graphical representation of {x ^{y)), integrated in the region y to y + 5 for different pion DAs. 
Left panel: Solid blue line — asymptotic; dashed green line — BMS [11]. Right panel: Dashed-dotted green 
line— CZ [3]; dashed-dotted-dotted red line— fiat-top DA [Eq. (B.l) with a = 0.1]. 



A. "Integral" sum rules 

We introduce now the integral derivatives in order to discuss the slope of the pion DA. 
We define these quantities in terms of a set of operators D given by 

X 

[dW^](x) = ^'(x), [D^'^^]{x) = ^{x)/x, [D^^)^]{x) = - l^dy. (3.2) 

J y 



Then, assuming (p{0) = 0, we get the iterative formula 

X 

[D(^^'^y,]{x) = ^ jdy[D^-^^]{y) . (3.3) 



Thus, each higher derivative within the set of the differential operators D^''^\ is stronger 
averaged with respect to x than the previous one. The usefulness of these derivatives derives 
from the fact that they can be applied to QCD sum rules which in general contain on their 
RHS singular contributions. In Appendix D, we elaborate on D^"^^ so that here we can focus 
our attention on the main properties of these derivatives. First, it is obvious that D^"^^ acts 
on a linear function as a differentiation operator, i.e., D^'^^ax = a. Second, assuming that 
the Taylor expansion of (p{x) at x = exists, one finds from (3.5) 

[D^^^\]{x) = ^'(0) + ^"(0)^ + O (l^^^^)) , (3.4) 

which is valid for any real u, as we explain in Appendix D. From the above equation, one can 
see that the defined operator D^'^^ reproduces at small x and/or large u the derivative of (p{x) 
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at the origin x = 0. Strictly speaking, using (3.5) and (D.4), one obtains \im[D^'^~^'^^ (p]{x) = 

</9'(0) (at fixed z/ G M) and lim [D^'^^'^^ ip](x) = (/''(O) (at fixed x). For this reason, we 

may appeal to Eq. (3.3) and call the variation range of the operator the "integral 

derivative" of (p. Having defined this operator, we can derive the following expression 



X 

1 



+ ^ _ j^^y)f^y^y^^)dy , (3.5) 







f{y^u^x) = ^^\ L^-\ (3.6) 



where 



T{v + l)y\ y 

for any real v (see Appendix D). As it is seen from Eq. (3.5), the function /(?/, v^x) acts as 
a "smooth projector" onto the vicinity of the origin of ?/, as one can appreciate from Fig. 10 
in Appendix D. 

By applying the operator [Z^^'^+^^j on both sides of the QCD SR given by (2.4), we obtain 
a new SR, viz.. 



(3.7) 



2\ 



[Z^^'+'Vpcrt] (x) e"^/*' rfs + [D('^+2)A<I>g](x,M2) + [D^''+^^A^y]{x,M 

+ [D(''+2)A$t](x, M^) + [D("+2)A$s](x, M^) . 

In order to achieve a better stability of this SR, we take into account an effective Ai-meson 
state that embodies the vr' and the Oi mesons and has the decay constant = 0.227 GeV 
and the mass Tn\^ = 1.616 GeV^. For the pion-decay constant and the continuum threshold, 
we use /,r = 0.137 GeV and Sq^° = 2.25 GeV^, respectively. These values were derived 
before from the corresponding two-point QCD SRs with nonlocal condensates, see Ref. [11]. 
However, because in our SR — cf. (3.7) — we have to resort to the LO expression for the 
spectral density (recall what we said about this issue before), we adopt a somewhat larger 
value of the threshold parameter: sq ~ s^^*^(l + as/n) = 2.61 GeV^. This is done for both 
the integral as well as the differential SRs, the reason being that we want to preserve the 
correct normalization of the pion DA.^ 

As usual, we study the SRs in the fiducial interval of the Borel parameter G 
[M^jjj, M^g^j^], where both terms, the continuum contribution and the nonperturbative one, 
each contributes about 1/3 to the whole SR (2.4). This induces an uncertainty of the order 
of (1/3)^ — )■ 10%. Moreover, the quantities to be calculated with the SR (the integral deriva- 
tives) should not (crucially) depend on the Borel parameter. Therefore, we should take 
care that this dependence is minimized. To achieve this goal, we attempt to minimize the 
root-mean-square deviation by varying the [D^'^~^'^'^lpai\{x) contribution in the fiducial Borel 



^ We thank A. P. Bakulcv for useful remarks on this point. 
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interval. On this account, we can average the M^-dependence of the pion DA contribution 
[first term in Eq. (3.7))] in order to get a more rehable form of the SR. Therefore, we write 

[D(^+2)y,f](x) = ([D('^+2)(^.](x,M2)) 

In Fig. 2 we show the M^-dependence of [D'^'^'^lPtj]{x, M"^), obtained from SR (3.7), for different 
values of [D^^'^LfAilix) and evaluating it for x = 0.5. The average value in the fiducial interval 
is [L'(3)^SR](^Q_5j ^ 4.8 ±0.5. 

6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 

0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 2: Dependence of [I?'^'^V7r](a;, Af^) on the auxiliary Borel parameter AI^ for x ~ 0.5. The solid blue 
line corresponds to [D^^^AiKa;) = 6.4, whereas the dashed lines refer to Lp a\\{x) = 6.7 (upper curve) 
and [D^^'^ipAilix) = 7.1 (lower curve), using in all cases x = 0.5. The two vertical lines delimit the fidelity 
region M^G [M2;^,M^,J. 

We discuss now the integral sum rules and their applications. The main contribution 
from the singularities in the SR for [D^'^^'^^lPt^]{x, M^) stems from the x-region around A = 
A^/(2M^). This is because we used a delta-ansatz model [cf. (A.l), (A. 2)] for the condensates, 
which implies that the nonperturbative contributions have delta-function and Heaviside- 
function behaved terms (A.3)-(A.8). Therefore, in order to take into account all NLC 
contributions, we should analyze the region x > 0.4 and use M^j^ > 0.6 GeV^ which 
corresponds to A < 1/3 < x. Moreover, the image of the operator D^'^^'^^ for z/ > 4 
is numerically very close to the result obtained with the differentiation method (see next 
subsection) — for any x. Thus, the integral SR (3.7) becomes close to the differential SR 
which we will consider in the next section. For these reasons, we analyze the constructed SR 
(3.7) for u = 0,1, 2, 3, 4 and x > 0.4 and present the results in Fig. 3 by the solid line that is 
inside the light gray strip bounded by the short-dashed lines. For the sake of comparison, the 
predictions for the asymptotic DA (dashed-dotted line) and the BMS DA bunch — obtained 
in the NLC SR analysis of Ref. [11] — (shaded band limited by long-dashed lines) are also 
shown. From this figure we see that our SR estimates for [D^''+^^ ^l^]{x) agree fairly well with 
the BMS model — see also Table I. This table shows estimates for the third-order integral 
derivative of the pion DA for x = 0.5, using (i) the sum rule given by Eq. (3.7) and (ii) the 
pion DA models we discussed above, and also flat-type DAs which we consider below. 

We turn now our attention to flat-type DAs. First, we compare the QCD SR result, 
obtained in (3.8), with what one finds with the flat-top model f^^^i^{x) ~ (x(l — x))" given 





[Z)(3)^^](0.5,M2) 

^ ^ 






^ - 






mm 


mfix 
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FIG. 3: We show the x-dependence of [D^''+'^"i ip^] {x) for the BMS bunch of pion DAs [11] (shaded green band 
within long-dashed hnes) in comparison with the SR resuh (3.8) (narrow gray strip) in ah three panels. The 
left panel shows the predictions for v — whereas those for 1^ = 1 and v ^ A arc shown in the middle and the 
right panel, respectively. The dashed-dotted line represents the asymptotic result [£)('^+-^'(p^'^](a;) = 6— 3a;/2''. 



in Appendix B. For this model, one has 

[^^'"^'V?Ri)](2;) > [I^("+'V'1(a;) > [D^''^^^^f']{x) for a < 0.1 

for any real v G +M and < x < 1 For the value a = 0.1, we find [D^^VfB*!)^^-^) = 227, 
which is much larger and far outside the range of values extracted from our SR. 

Second, we consider a particular flat- type pion DA which is provided by the AdS/QCD 
correspondence in the holographic approach — see, for instance, Refs. [40-43]. In that case 
one has a = 0.5 yielding [D^'^^if'^°^]{0.5) = 14. 

Third, we study an alternative fiat-like pion DA which results from the Gegenbauer 
expansion of unity by retaining only the first few harmonics. One obtains 

3 

^f^Ux) = Qxxycli\2x - 1) , ^(^^ + 3) (3 9) 

Y-(3.9)V ; ^ 2n\ ^3 2n + l 2n + 2 ^ ' 

with a profile shown in Fig. 4 in comparison with the models already mentioned: BMS — solid 
line; CZ — long-dashed blue line; fiat-top DA given by Eq. (B.l) — dotted red line; fiat-like 
DA given by Eq. (3.9) — dashed-dotted green line. 




FIG. 4: Comparison of selected pion DA models. Solid blue line — central line of the BMS bunch [11]; 
long-dashed blue line — CZ model[3]; dashed-dotted green line — flat-like model given by (3.9), dotted red 
line — flat-top DA from Eq. (B.l) with a = 0.1. All DAs are normalized at the same scale ~ 1 GeV^ with 
or without evolution [11]. 
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It is interesting to notice that using the pion DA given by expression (3.9) to calculate the 
pion-photon transition form factor according to Radyushkin's expression (24) in Ref. [33], 
one actually reproduces the gross features of his results. This proves that the inclusion of the 
first few Gegenbauer polynomials in (3.9) does not affect the result obtained with (p{x) = 1 
in the momentum range < 40 GeV^ in a crucial way. This type of DA, i.e., (3.9), yields 
for the integral derivative the value [D^^^Lp^^^g-^](0.5) = 22.5, which is much larger than the 
range of values determined via our SR. This holds true also for the other two fiat-type 
DAs considered above. Recalling that the leading-order QCD sum rules with the minimal 
Gaussian model for the nonlocal condensates provide much smaller values of the integral 
derivative of the pion DA, one may conclude that it is very difficult to reconcile fiat-type 
pion DAs with SR (3.7). 

On the other hand, also the CZ model yields third-order integral derivatives which are 
incompatible with the values derived from our SR (3.7) — see Table I. Note that a similar 
statement also applies to the pion DA proposed in [44], which employs a Brodslcy-Huang- 
Lepage ansatz for the k^-dependence of the pion wave function — see Table I. The upshot of 
this table is that the SR for the integral derivative of the pion DA is fulfilled by the BMS 
bunch, whereas fiat-type DAs have no overlap with the estimated range of values. The same 
is true for the CZ model. 



TABLE I: Results for the third-order integral derivative for x = 0.5 and the (usual) derivative of 
the pion DA, using different SR approaches (first three rows) and pion DA models (six last rows). 





Approach /Model 


"Integral" derivative [D''^'>ip^]{0.5) 


Derivative (/?^(0) 


1 


Integral SR (3.7) 


4.7 ±0.5 


5.5 ± 1.5 


2 


Differential SR (3.14) 




5.3 ±0.5 


3 


SR (3.14) with smooth NLC (3.18) 




7.0 ±0.7 


4 


BMS bunch [11] 


5.7 ± 1.0 


1.7 ±5.3 


5 


Asymptotic DA 


5.25 


6 


6 


CZ DA [45] 


15.1 


26.2 


7 


DA from [44] 


14 





8 


Flat-type DA, Eq. (3.9) 


22.5 


72 


9 


flat-top DA (Eq. (B.l), a = 0.1) 


227 


>6 



We close this subsection by considering the usual derivative v^^(O) of the pion DA which 
encapsulates the key characteristics of the pion DA at small x. This quantity can be extracted 
from Fig. 3, where we have plotted the results for [D^'^^ (f^^]{x) we obtained with the SR (3.7) 
for different values of u. To determine </?^(0) one can use expansion (3.4) 

[D^^^'^^f']{x) ^ ^'M + ^:(0)^ (3.10) 

and subtract the second derivative for which the asymptotic value </?"(0) = —12(6) is used. 
The involved error ±6 was estimated by varying (3.10) within the narrow grey strip in Fig. 
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3. To be more specific, one obtains 

^:(0) ^ [D(^+^Vf ](x) - ^':(0)^ ^ [^('^^^Vf ](x) + I = 5.5 ± 1.5 (3.11) 

for any 0.4 < x and < u < 4. The error in (3.11) is a combination of tlie uncertainties 
originating from SR (3.7) and tfie error in tlie determination of ^p"{0). 

The above finding can be compared with what one obtains for the BMS and the CZ model 
(displayed in Table I) 

^1(0) = 6 [1 + 6 + 15 a^'] = l-07t'^m = "3.61 ^ 6.95 (3.12) 

and 

(^UO) = 6 [l4-6aP] -26.2 , (3.13) 

respectively, using in both cases the normalization scale /x^ ~ 1 GeV^. Quite analogously 
to the integral derivative, the CZ model yields also for the standard derivative much larger 
values than those estimated in (3.11). By contrast, the pion DA model proposed in [44] — 
though it provides a similarly large integral derivative like the CZ DA — has a usual derivative 
at the origin which is zero due to the strong exponential suppression of this DA in the vicinity 
of the origin. 



B. DifTerential sum rules 

Another way to study the behavior of the pion DA in the small-x region is provided by 
the differentiation of the SR (2.4), which yields 

^UO, M') = (l - e-^o/*^^) + 18^5$' - fl . (3.14) 

We shall evaluate this SR for the threshold value Sq = 2.61 GeV^, recalling that we employing 
a LO expression for the spectral density. 

Using the simplest delta- ansatz model for the condensates (cf. (A. 2)), only one nonper- 
turbative term survives, namely, the four-quark-condensate 

^' = -^^A^s{x,M') , (3.15) 
18As dx x=o 

where A^si^, M'^) is represented by Eq. (A. 9) — see Fig. 8 in Appendix A. The vector-quark 
condensate (V) and the gluon condensate (G) — recall Eq. (2.4) — give zero contributions in 
the region a; < A, where A = A^/(2M^) > with = 0.4 GeV^. On the other hand, 
the antiquark-gluon-quark condensate (T) amounts to a vanishing contribution in the region 
X < min { A, 1 — 2A}. For > 0.4 GeV^ (or equivalently A < 1/2) this contribution can 
also be neglected. 

Even if we assume a behavior of the various condensates differing from the delta-ansatz 
model, using, for instance, a smooth model like (3.18) for the scalar quark condensate (which 
implies a decay at large distances not slower than the exponential decay — see further below), 
the (V), (G), and (T) terms give only a small non-zero contribution in the small-x region. 
Therefore, these terms can be neglected in first approximation in both mentioned models. 
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FIG. 5: The distribution of vacuum quark virtualities in the smooth model (3.18) for a particular choice of 
the intrinsic parameters as explained in the text. 

In this section, we shall study the differential SR (3.14) using the mentioned models for the 
scalar condensate. 

The calculation of the nonlocal version of the four-quark-condensate contribution involves 
only the scalar quark condensate. The four-quark contribution A$s(3;; was obtained in 
[4] and can be used with various types of vacuum-quark distributions /5(a). Accordingly, 
the probability for finding vacuum quarks with very large virtualities is very small, as one 
can see from Fig. 5. Therefore, the distribution function /5(a) ~ diminishes in the region 
a > M^, the latter region corresponding to the values > M^^^ > 0.6 GeV^ of the Borel 
parameter. On the basis of this result, one [46] can use the method of Ref. [4] to obtain 
the model-independent expression (A. 9). Then, one finds for the four-quark-condensate 
contribution to the SR (3.14) the following expression 

$' = Hda^^ = {qqy rz^Ms{z^)dz^ . (3.16) 
Jo a Jo 

We see from this equation that the nonperturbative contribution to the SR is mainly due to 
the scalar-quark condensate at large and moderate distances z"^ ~ 4/(/c^). 

For the concrete evaluation of the differential SR we again employ the same criteria as 
already used in the integral SR for both the continuum and the nonperturbative terms. 
Applying these criteria, the low boundary for the Borel parameter turns out to be very 
small, viz., M^^^ = (0.3 — 0.4) GeV^. Though this low value is consistent with the applied 
criteria, it is not in good agreement with standard QCD sum-rule approaches in which a 
higher minimum Borel-parameter value is used. Therefore, we use M^^^ = 0.6 GeV^, a value 
also employed in the QCD SR for the pion DA in Ref. [11] in connection with the moments 
of the pion DA. Keep also in mind that using a lower value of the Borel parameter would 
cause the decrease of the first derivative of the pion DA at the endpoints. To continue, we 
define the derivative v^^(O) = (</?^(0, M^)) in terms of the mean value in the fiducial interval 
G [M^jjj, M^^^] following the definition on the RHS of (3.8). This helps minimizing the 
sensitivity of 9?^(0) on the choice of the Borel parameter by means of the variation of the Ai 
contribution (p'^_^{0). 

The delta-ansatz /5(a) = 6{a — Xg/2) leads to a simple expression for the nonperturbative 
contribution to SR (3.14), notably, 

= -^delta = 4 • (3.17) 
Q 

Note that we defined the Ai-meson contribution </?^^(0) = 6.7 by means of the minimum of 
the root mean square deviation. On the other hand, the dependence of ip'^{0,M'^) on the 
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FIG. 6: Left panel: The solid curve shows the M^-dependence of (p'^{0,M'^) in the SR with the Ai- 
meson contribution (p^^(O) = 6.7. The broken lines represent (^^(0,M^) for <(5^^(0) — 6.1 (lower line) and 
fAii^) — '''■^ (upper line). Right panel: The solid curve illustrates the Af^-dependence of (^^(0, M^) in the 
SR which employs the smooth quark condensate model with a Ai-meson contribution given by tp'j^_^ (0) ~ 6.8. 
The broken lines denote (p'^{0,A'P) at Lp'j^_^{0) — 6.1 (lower line) and (p'j^_^{0) — 7.5 (upper line). 



Borel parameter for the delta-ansatz model is controlled by Eq. (3.14) and is shown in 
the left panel of Fig. 6. Thus, the average value of the pion DA derivative in the fiducial 
Borel interval is v^^(O) = 5.3(5) — see Table I. 

We go forward and discuss the consequences of the smooth model for the quark-virtuality 
distribution in the differential SR. 

Though the delta-ansatz model is useful, because of its simplicity, there is an indica- 
tion from the heavy-quark effective theory [47] that in reality the quark-virtuality distribu- 
tion fs should be parameterized in a different way as to ensure that the scalar condensate 
decreases exponentially at large distances. Moreover, in order that the vacuum matrix ele- 
ment (g(D^)^g) exists, the quark-virtuahty distribution /s(a) should decrease faster than 
any power l/a^~^^ as a — > oo [4]. For this reason, the authors of [36, 48] suggested a two-tier 
model for fs which has a smooth dependence on the quark virtuality a, namely, 

where Kn{z) is the modified Bessel function. This model, the so-called "smooth model", 
depends on two parameters A and a that parameterize, respectively, the long- and short- 
distance behavior of the nonlocal condensates [36]. For large distances | z\ = \J —z^ this 
model leads to the asymptotic form 



I 9(2n-l)/2 

VA A„(2Acrj 

It is instructive to consider a purely exponential decay of the quark-virtuality distribution 
and study its influence on the quark condensate. This can be realized in the model of [36, 48] 
by choosing n = 1, whereas the second parameter A = 0.45 GeV can be taken from the 
QCD SRs for the heavy-light meson in heavy quark effective theory — see [47, 49]. The two 
parameters n and A are responsible for the large- 2; behavior of the scalar-quark condensate, 
cf. Eq. (3.19). The third parameter cr^ = 10 GeV~^ is defined in terms of the parameters 
n. A, and via the following equation: 

" f r A \A A K„+i(2Aa) \\ 

a fs{a;A,n,a) da = ———— = —, (3.20) 

a Kn[ZAaj 2 
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FIG. 7: Dependence of the derivative <y5^(0) on the parameter n of the smooth scalar condensate model (3.18) 
for A = 0.3 GeV (dashed-dottcd red line), A = 0.45 GeV (solid blue line), and A = 1 GeV (dashed green 
line). Here, we show only the central value of the derivative V'^(O) we obtained from the SR analysis (3.14). 
The black dot symbol marks the position which corresponds to the model parameters n = 1, A = 0.45 GeV, 
and 'y3^(0) = 7.0(7), that corresponds to the third row in Table I. 



which we evaluate for the value of the nonlocality parameter = 0.4 GeV^. The main effect 
of using a smooth model for the quark- virtuality distribution relative to the Gaussian form, 
fs{ci) = 6{a — A^/2), is the induced increase of the nonperturbative contribution to the SR, 
so that 

^smooth = I da = > ^<^e,ta • (3-21) 

We analyzed the SR (3.14) for this model using a particular choice of its parameters, 
notably, fs{a;A = 0.45 GeV, tt, = 1,0"^ = 10 GeV~^) and determined the dependence of 
ip'^{0,M^) on the Borel parameter M^. The result is shown graphically in the right panel 
of Fig. 6. The average value of the derivative ip'^{0,M'^) in the fiducial Borel interval is 
(p'j,{0) = 7.0(7). Thus, the nonperturbative contribution $gi„ooth' obtained from the smooth 
model, is approximately two times larger Xthan the analogous contribution ^d^ita 
delta ansatz: ^smooth ~ 2.3 $dcita- addition to this result, marked by a black dot, we show 
in Fig. 7 the dependence of V'^(O) on the choice parameters n and A of the smooth model. 
From this picture and the relation (3.16), we may come to the conclusion that, choosing a 
model for the condensate that has a slower decay at large distances (small n or A), may 
cause an increase of the nonperturbative contribution to the SR (2.4) and entail also an 
increase of the value '/'^(O). On the other hand, choosing a model for the condensate with a 
faster decay at large distances (large n or A), may lead to a decrease of the nonperturbative 
contribution to the SR (2.4) and therefore to a decrease of the value (p'j,{0). 

In the last column of Table I we collect the values of the (usual) pion DA derivative at 
X ~ 0, using different SR approaches (first three rows) and pion DA models (last six rows). 



IV. CONCLUSIONS 



The current investigation was partly motivated by the recent results obtained by the 
BaBar Collaboration [30] on the pion-photon transition form factor which indicate an unex- 
pected growth of this quantity with above ~ 10 GeV^ up to the highest momentum value 
of 40 GeV^ measured. As it was pointed out in [27] (see also [31, 32]), such behavior is 
impossible for endpoint-suppressed pion DAs and brings into play a fiat-type profile for the 
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pion DA, as proposed in Refs. [33, 34]. Therefore, it appears to be of crucial importance to 
have a theoretical tool in our hands able to reveal the particular characteristics of pion DAs 
precisely in the kinematic endpoint region. 

The method mostly used in the past to extract such information employs the inverse 
moment — extensively discussed in [11, 22, 23, 25, 37]. In the present work we pro- 

posed another, more direct way, to access the endpoint characteristics of the pion DA which 
makes use of its derivatives. The first step was to define the notion of "integral" derivatives 
[D^'^^ (pTr]{x) by means of an appropriate operator D^'^^ [cf. (3.5)]. Next, we formulated an 
"integral" sum rule, Eq. (3.7), for these quantities by taking into account an effective Ai- 
meson contribution. Using this sum rule, we determined the range of values of the integral 
derivatives and displayed it in the first entry of the second column of Table I showing it also 
graphically in Fig. 3. We also calculated the integral derivative for different characteristic 
pion DA models and listed its value in the same Table. The list of pion DAs includes the 
BMS, the CZ, the asymptotic, and two options for fiat-type DA models, one given by Eq. 
(B.l) with a = 0.1 — "fiat-top" model — the other being a fiat-like model parameterized in 
terms of Eq. (3.9). We also analyzed the usual derivative of the pion DA in the vicinity of 
the origin. This procedure is helpful in revealing the fine details of the pion distribution 
amplitude in the endpoint region — see Table I for the results. 

Our findings can be summarized as follows. First, at the end of Subsec. Ill A, we applied 
the result for the integral derivatives [D^'^~^'^'>Lp.,^]{x) in order to reproduce the usual derivative 
of (fTrix) at the origin x = 0. The result v^^(O) = 5.5 ± 1.5 is shown in the third column of 
Table I in comparison with some characteristic models for the pion DA. Second, in Subsec. 
(IIIB), we studied the derivative of the pion DA v^^(O) in terms of the differentiation of the 
SR (2.4) that leads to the differential SR (3.14). It turns out that the only nonperturbative 
content in this SR is mainly defined by the scalar quark condensate. The nonperturbative 
contribution is proportional to the second inverse moment (3.16) of the distribution fs{o:) 
of the vacuum-quark virtualities and is defined by the behavior of the quark condensate at 
large and moderate distances between the vacuum quarks. The results for the derivative 
of the pion DA v^^(O) are shown in Table I for the delta-ansatz model (A.l) and also for 
the two-tier smooth model (3.18) with the model parameters n = 1, A = 0.45 GeV. The 
dependence of the derivative of the pion DA v^^(O) on the choice of the model parameters n 
and A is illustrated in Fig. 7. 

As we see from Table I, the results from the differential (3.14) and the integral SRs (3.7) 
agree with each other. The integral derivative of the pion DA, based on a new SR derived 
in this work, remains smaller than the asymptotic value and overlaps with the range of 
values determined with the BMS bunch of pion DAs, while there is no agreement with the 
CZ DA and the fiat-type models considered. The same conclusions can be drawn also for 
the usual derivative of the pion DA, which follows from the differential SR (3.14). It is 
worth mentioning that employing the integral and the differential sum rules, we found that 
the leading-order QCD sum rules (2.4), which employ the minimal Gaussian model for the 
nonlocal condensates, cannot be satisfied by flat-type pion distribution amplitudes. 

Given that an increasing behavior of the scaled pion-photon transition form factor can 
only be achieved with fiat-type pion DAs, the independent experimental confirmation of 
this effect, e.g., by the BELLE Collaboration, becomes extremely crucial for our theoretical 
understanding of basic QCD exclusive processes. 
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Appendix A: Expressions for the nonlocal contributions to the sum rules 



In Sections II and III, we used the following expressions for the vacuum distribution 
functions 



fs{u) = S{u-Xy2) ■ fy{u) = 6'{u-Xl/2) 
/t,(«i, «2, «3) = S{ai- A5/2) 6 (^2 - Aj/2) 6 («3 - Aj/2) . 



(A.l) 
(A.2) 



The meaning of these expressions and their connection to the initial NLCs has been discussed 
in detail in Refs. [4, 7]. The contributions to the QCD SR (2.4) with nonlocal condensates, 
A$r(x, M^), associated with these expressions, are shown below. Here, and in what follows, 
we use A = Xl/{2M'^), A = 1 - A. 
Then we obtain 



A<l>5 (x, M^) 

A^v (x,M2) 
A$T {x, M^) 



A$T2 (x, M^) 



A<^T, {x, M^) 



A$G {x, M^) 



As 18 
AA2 



{9 {x > A > x) x[x + {A - x) In (x)] + (x x) + 



1 > A)^ (A > X > A) [A + (A - 2xx) ln(A)] } , (A.3) 

(A.4) 

A<I>T, (x, M^) + A<I>T, (x, M^) + A$T3 (x, M^) , 



^(x5'(x-A) + (x^x)) 



' M4 



[6{x - 2A) - 6{x - A)] 



x9{x > A)9{x > 3A - i; 



X 

A 



3x 
A^ 



6 



^(1 > 2A) + e(2A > X 
1 + X 



A 



+ (x — !■ x) 



(A5) 



4As__ | ^(x > 2A) - ^(2A > x)^(x > A)^(x > 3A - 1) 

1 2x - 4A 



X -7- X 



3Asx r 



M^AA i 
-|- (x — )■ X 
(asGG) 



AA2 

^(2A > x)^(x > A)^(x > 3A - 1) 



9 _ 1 _ ^' 
A A 



247rM2 



(5(x- A) + (x ->x)) . 



(A.6) 

(A.7) 
(A8) 
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In the above equations, we used the abbreviation As = ^^f^^ (gg)^, while for the quark and 

81 

the gluon condensates the standard estimates [50] as{qq)'^ = 1-83 ■ lO"'' GeV^, — = 

2{qq) 

been adopted. In Sec. Ill B, we applied not only the delta ansatz, discussed above, but also 
the smooth model, proposed in Refs. [48, 51], which is based on an exponential decay of 
the condensate (3.19). For this reason, we use the model-independent expression for the 
four-quark contribution [4, 46] to obtain Eq. (3.16) 



0.0012 GeV^ and = = 0.4 GeV^, normalized at /i^ ^ 1 GeV^, have 



oo 



A<S>s{x.M^) = ^ ljda,da,Js(a,)fs(a,) (A.9) 



00 

x6' (Ai - x) 
+ (x -> x) , 



where Aj = ai/M"^, Aj = 1 — Aj, and x = 1 — x. 




FIG. 8: Typical diagram for the four-quark nonlocal-condensate A$s(a;,Af^) encoding nonperturbative 
input in (2.4). 



Appendix B: Modeling the pion DA 

In our work we compared the results following from the differential SR (3.14) for the 
usual derivative of the pion DA, v^^(O), with those obtained with the help of the integral 
derivative, [D'^'^'^lPt^]{x)^ i.e., (3.7), evaluating them for the BMS DA bunch, the CZ model, 
and a couple of flat-type pion DAs. 

The BMS bunch of pion DAs was obtained [11] as the sum of the first three terms in 
the Gegenbauer-polynomial expansion (1.3). Their Gegenbauer coefficients were obtained 
from an analysis of the SR (2.4) for the (^^^)-moments (A^ = 0, 1, 2, 3, 4, 5) with the nonlocal 
condensate contributions being presented in Appendix A. It was found that only the two first 
Gegenbauer coefficients 02 and 04 give sizeable contributions to the (^^^)-moments, whereas 
the higher Gegenbauer terms give merely tiny contributions. For this reason, the BMS bunch 
DAs are two-parameter models. The central values of a2 and 04 for the whole BMS bunch 
in the 02,04 space are a^^^ = 0.187 and 04^^ = —0.129 using as a normalization scale 
/i^ - 1 GeV^ [11]. 

On the other hand, the CZ pion DA contains only the first nontrivial Gegenbauer poly- 
nomial corresponding to the coefficient 02 [3]. It was derived from QCD SRs for the (^^^)- 
moments {N = 0, 1,2), using local condensates. For the sake of consistency, we use here a 
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value of = 0.56 obtained after evolution to the normalization scale //^ ^ 1 GeV^, see for 
more details [22]. 

One of the two flat-type models considered in this paper is defined by 



For this model, we find the following expressions: 



These expressions can be generalized to any real differentiation index u G +IR for < x < 1 
to get 

[^^''^'VcrdK^) » [D^-'^'^^nix) > [Z^(^+^Vf ](^) for « < 0.1 , 
227 > 5.25 > 4.7(5) 

for the particular values a = 0.1, u = 1 and x = 0.5 , 

[^^^^'V?B^)](^) > [D'^'^'^vni^) > [D(^+^Vf ](^) for « < 1 , 

14 > 5.25 > 4.7(5) (B.2) 
for the particular values a = 0.5, u = 1 and x = 0.5 . 

The other fiat-type model, given by Eq. (3.9), was already discussed in the text. 
Appendix C: The spectral density 

The spectral density Ppci-t{x) with NLO accuracy was obtained in [4, 7] and was found to 



be 

1 



Ppert{x) — 3XX 



l + ^CFi5-- + \n\x/x) 



27r2 



(C.l) 



For our numerical calculations, which take into account radiative corrections, we use the 
following function 

The x-dependence of the function A(z/, x) is illustrated in Fig. 9 for different values of 
solid blue line — u = 0, dashed green line — u = 1, and dashed-dotted red line — u = 4. 

Appendix D: Properties of the integral operator 

In order to ensure a weak dependence of the results on the particular model for the con- 
densates adopted, and in order to include all condensate contributions, one has to construct 
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the SR by integrating the pion DA SR (2.4) over a large enough interval of x. For this rea- 
son, we introduced in Sec. Ill A the integral derivatives [D^^Vll^); with the two lowest-order 
ones being given in terms of Eq. (3.2). The next higher derivative reads 

t 

Assuming Lp{Q) = 0, we find 

X 

[D("+iV](a;) = lj'^y [^^"V](a;) • (D-2) 



To obtain an expression for [D^'^^'^^ (p](x) for any n > and n G N, one has to rearrange the 
integration order to get 

\t / 

This can be readily generalized to the expression (3.5) using (3.6) that allows one to establish 
this transformation for any real i^, z/ G M. The dependence on y of the function f{y,u,x) 
is shown in Fig. 10 for x = 0.6 and u = 0,1,4. As one sees from Eq. (3.5), the function 
f{y, u, x) acts as a smooth projector to the area around the origin of y. 

If a Taylor expansion exists for (p{x) at x = 0, then by applying (3.5) one finds 

oo „ 

p-vi(x) = ^'(0) + ^"(0)^ + j: ^'""'(0) ;^ . (D.4) 

71=2 ^ ' ^ ' 

Using (3.5) in combination with (D.4), one can obtain the following properties of the operator 
lim[L)(^V](a;) = hm [D^^VK^;) = ^^'(0) and D^'^) {ax + hx^) = a + l^-^'hx. Moreover, 

x— >0 I/— >oo 

from (D.4) it follows that the introduced operator \D'^^^i^{x) reproduces at small x and large 
V the derivative of '^{x) at a; = 0. 
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